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[57] ABSTRACT 

Young’s fringes are produced from a double exposure 
image of particles in a flowing fluid by passing laser 
light through the film and projecting the light onto a 
screen. A video camera receives the image from the 
screen and controls a spatial light modulator. The spa- 
tial modulator has a two dimensional array of cells the 
transmissiveness of which are controlled in relation to 
the brightness of the corresponding pixel of the video 
camera image of the screen. A collimated beam of laser 
light is passed through the spatial light modulator to 
produce a diffraction pattern which is focused onto 
another video camera, with the output of the camera 
being digitized and provided to a microcomputer. The 
diffraction pattern formed when the laser light is passed 
through the spatial light modulator and is focused to a 
point corresponds to the two dimensional Fourier trans- 
form of the Young’s fringe pattern projected onto the 
screen. The data obtained from the second video cam- 
era may be processed by the computer to find the corre- 
lation peak locations, areas, and other data of interest at 
the portion of the film through which the laser beam 
was projected. The film may be indexed such that the 
beam is passed through the film to cover the entire area 
of the double exposure image recorded on the film and 
thereby provide information concerning particle dis- 
placement and direction of the displacement at each 
position in the film. 

24 Claims, 7 Drawing Sheets 
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OPTICAL CORRELATOR METHOD AND 
APPARATUS FOR PARTICLE IMAGE 
VELOCIMETRY PROCESSING 

5 

This invention was made with U.S. Government 
support awarded by the Department of the Army 
(DOD) and NASA grand number(s): DOD 

#DAAL03-86-K0174 and NASA #NAG3-718. The 
U.S. Government has certain rights in this invention. 10 

A portion of the disclosure of this patent document 
contains material which is subject to copyright protec- 
tion. The copyright owner has no objection to the fac- 
simile reproduction of the patent document or the pa- 
tent disclosure, as it appears in the Patent and Trade- 15 
mark Office patent file or records, but otherwise re- 
serves all copyright rights whatsoever. 

FIELD OF THE INVENTION 

This invention pertains generally to the field of fluid 20 
velocity measurements and particularly to particle 
image velocimetry techniques. 

BACKGROUND OF THE INVENTION 

Significant recent interest in full-field velocity imag- 25 
ing techniques has focused on speckle velocimetry or 
particle image velocimetry techniques. See, e.g., Ro- 
land Meynart “Instantaneous Velocity Field Measure- 
ments in Unsteady Gas Flow by Speckle Velocimetry,” 
Applied Optics, Vol. 22, No. 4, 15 Feb. 1983, pp. 30 
535-540; Ronald J. Adrian, et al., “Development of 
Pulse Laser Velocimetry (PLV) for Measurement of 
Turbulent Flow,” Symposium on Turbulence, Univ. of 
Missouri, Rolla, 1984, pp. 170-186; C. C. Landreth, et 
al., “Double Pulse Particle Image Velocimeter With 35 
Directional Resolution for Complex Flows,” Experi- 
ments in Fluids, Vol. 6, 1988, pp. 119-128. Despite the 
different terminology of some authors, most of these 
techniques are similar. The techniques evaluate particle 
velocities from double exposure photographs separated 40 
by a known time differential. For conditions in which 
the particle density is large, a common analysis tech- 
nique for the double exposure particle photographs is to 
pass an unexpanded laser beam through the double 
exposure image. The particles which fall within the 45 
beam scatter the beam, with particle image pairs pro- 
ducing the well known pattern of Young’s fringes. De- 
pending on the conditions under which the particle 
images were produced, the sets of particle pairs within 
the probe beam may exhibit a range of apparent dis- 50 
placements. Average displacement and direction infor- 
mation may be obtained by analyzing the Young’s 
fringe pattern for the fringe spacing and fringe angle. 
Information concerning the details of the variance of 
the displacement of the velocity pairs within the probe 55 
beam from the mean displacement is carried by the 
modulated shape of the Young’s fringes, but is difficult 
to extract directly. See, K. Hinsch, et al., “Fringe Visi- 
bility in Speckle Velocimetry and the Analysis of Ran- 
dom Flow Components,” Applied Optics, Vol. 23, No. 60 
24, Dec. 15, 1984, pp. 4460-4462. 

An alternative processing technique has been pro- 
posed which acquires the two dimensional spatial Fou- 
rier transform of the Young’s fringe pattern, providing 
correlation peaks in the frequency plane at the fre- 65 
quency and angle corresponding to the fringe pattern of 
Young’s fringes. See, Ronald J. Adrian, et al., supra. 

The shape of the correlation peak, as well as the overall 
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background level, can be calculated to provide quanti- 
tative estimates of the velocity variance at the particular 
laser beam probe point. 

A typical arrangement for carrying out the forgoing 
technique includes a laser beam positioner, the image of 
interest, and some type of high resolution image capture 
device, such as a 512 x 512 charge coupled device 
(CCD) array. For each beam position corresponding to 
a position in the flow field, a fringe pattern is produced, 
recorded, and numerically transformed. However, nu- 
merical transformation of 512 x 512 arrays can be time 
consuming, even for relatively specialized processors. 
A typical field may require 10,000 (100 by 100) interro- 
gation points, and an unsteady flow may require many 
images, depending on the levels of detail required. Con- 
sequently, a time consuming series of digital transforms 
are required. 

SUMMARY OF THE INVENTION 

The present invention carries out the two dimen- 
sional spatial Fourier transforming of the Young’s 
fringes utilizing optical techniques. Relatively inexpen- 
sive components are utilized, and the required digital 
processing time can be substantially reduced from that 
required by present processing techniques. 

In carrying out the present invention, a double expo- 
sure image is captured, as on a photograph, of particles 
in the flowing fluid in a conventional manner. The pho- 
tograph is then projected, using laser light, on a position 
by position basis onto a screen and the image on the 
screen is viewed by an image capture device, such as a 
video camera. The two dimensional image of the fringes 
is converted to an electrical signal by the video camera 
and is used to drive a two dimensional transmissive 
spatial light modulator such as a liquid crystal array. 
Collimated light from a second laser beam is then passed 
through the optical modulator, and utilizing proper 
optical components, the light from the modulator is 
focused onto a viewing screen or another video camera. 
The far field pattern which is focused on the camera is 
analogous to the result of a digital two dimensional 
Fourier transformation of the original fringe pattern. 
When the far field correlation pattern is recorded by the 
video camera, the signal from the video camera may be 
digitized and supplied to a computer for processing to 
find the correlation peak locations, areas, and other data 
of interest. This is repeated for each point on the film 
until the entire film image has been covered. In this 
manner, a two dimensional Fourier transform of the 
fringe pattern is provided at rates which are at least 
equal to or faster than specialized digital processors and 
at a significantly lower cost. 

Further objects, features, and advantages of the in- 
vention will be apparent from the following detailed 
description when taken in conjunction with the accom- 
panying drawings. 

BRIEF DESCRIPTION OF THE DRAWINGS 

In the drawings: 

FIG. 1 is a schematic side view of apparatus for ob- 
taining double exposure photographic images of parti- 
cles flowing through a test cell. 

FIG. 2 is a top view of the apparatus of FIG. 1. 

FIG. 3 is a schematic view of the apparatus of the 
present invention for carrying out the optical two di- 
mensional Fourier transformation of the Young’s fringe 
patterns and carrying out the analysis thereof. 
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FIG. 4 shows the typical output of a processor for a 
phantom fringe pattern, showing a central “DC” spot 
and off axis correlation peaks corresponding to the 
input fringe angle and displacement. 

FIGS. 5-7 shows three different fringe orientations 5 
and the respective convolutions. In each figure, the 
figure marked “A” is the input Young’s fringes as im- 
aged on a video monitor and provided to the optical 
modulator and “B” is the processor output after low 
frequency subtraction. 10 

FIGS. 8A and 8B are views as in FIGS. 5-7 but 
utilizing a lower frequency fringe spacing. 

FIG. 9 is a graph illustrating a comparison of a slice 
through the correlation peak of the simulated fringe 
pattern of FIG. 5A digitally transformed using a 256 x 15 
256 2-D FFT routine and a similar slice for the proces- 
sor output of the present invention. 

DETAILED DESCRIPTION OF THE 

INVENTION 20 

The present invention utilizes double exposure photo- 
graphic media image frames produced in a conventional 
manner to capture the information on particle displace- 
ment in the flowing fluid. FIGS. 1 and 2 illustrate the 
apparatus utilized to obtain these film exposures, with 25 
FIG. 1 being a side view and FIG. 2 being a top view. 
Light from a laser 11 (e.g., a 4 Watt argon-ion laser 
from Coherent Inc.) is passed through a shutter 12 . 
(e.g., a computer controlled acousto-optic modulator 
from Interaction Inc.) to cylindrical lenses 14 which 30 
produce a wide but very thin beam 15 of the monochro- 
matic light from the laser. The cylindrical lenses are 
selected for appropriate focal length and are available 
commercially (e.g., from Oriel Corp.) This beam is 
passed through a test cell 16 through which fluid is 35 
flowing. The test cell may be formed with a transparent 
wall or walls formed of a suitable material such as glass 
or acrylic plastic. The exiting beam is absorbed by a 
beam dump 18 in a conventional manner. The fluid 
flowing through the cell 16 has small particles therein 40 
which fluoresce or reflect the laser light. Exemplary 
particles are 1-2 micron titanium dioxide. A camera 20 
(e.g., a conventional 35 mm camera) is positioned to 
capture the image of the fluid flowing through the test 
cell 16 and faces the broad side of the beam 15. To 45 
obtain an exposure, the shutter of the camera 20 is 
opened and the shutter 12 is quickly opened and closed 
twice in succession, with a precise and known time 
differential increment between the times at which the 
shutter is opened. The particles in the fluid moving 50 
through the ceil 16 will be displaced a small distance in 
the fluid, evidencing the local motion of fluid during the 
time between openings of the shutter 12. Thus, each of 
the particles in the fluid will show up twice on the film 
exposed in the camera 20, separated by a small distance 55 
which is proportional to the velocity of the fluid at the 
particles, and with the direction of displacement of the 
particles indicating the direction of motion of the fluid 
with respect to the plane of the laser beam 15. 

The film from the camera 20 is then developed, 60 
contact printed onto film for a positive image, and 
placed in position in a holder 30 shown in FIG. 3. A 
stepper motor 31 (e.g., 0.00001 inch step displacement 
motors from Daedel, Inc.) drives the holder in two 
dimensional motion. Light from a laser 32 (e.g., 0.5 mW 65 
helium-neon laser from Uniphase, Inc.) is focused in a 
narrow beam 33 and passed through a portion of the 
film at each position of the film. The light passing 
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through the film is projected on a ground glass screen 
35 (e.g., Rolyn Optics ground glass). The image on the 
screen 35 is recorded by a video camera, 36 such as a 
vidicon camera (Xybion Corp.), which provides an 
electrical signal on its coaxial cable output lines 37 
which correspond to the image on the glass screen 35 in 
the conventional raster scan. The stepper motors (2 
axes), under control of a microcomputer system 40, are 
controlled to move the film on the holder 30 in step 
increments in a raster fashion so that the entire film has 
the beam 33 from the laser passed through it, with the 
resulting image at each position being recorded by the 
video camera 36. The microcomputer system 40 may 
consist of a microcomputer (e.g., a PC AT or clone 
such as Compuadd) and a frame grabber (e.g., Imaging 
Technologies FG-100), with a digital input/output 
board to drive the stepper motors. 

At each position on the film at which the laser beam 
passes through the film, the light projected from the 
film onto the ground glass screen 35 will form a set of 
Young’s fringes. The output signal on the line 37 from 
the video camera 36, which is the electrical signal corre- 
sponding to the two dimensional raster scanned image 
of these fringes, is utilized to control a two dimensional 
transmissive spatial light modulating device 45. This 
device may be an inexpensive modified liquid crystal 
television screen (e.g., Radio Shack 16-156 LCTV). 
The cells of the (e.g., 256 x 256) array of transmissive 
elements in the modulator 45 are controlled by the 
signal on the line 37 from the video camera 36 to attenu- 
ate the light passed through them in relation (e.g., in- 
versly proportional) to the brightness signal for the 
corresponding pixel of the image recorded by the video 
camera 36. The image presented on the modulator 45 is 
thus essentially the Young’s fringes which appear on the 
ground glass 35. 

A laser 47 (e.g., Siemens 5 mW helium-neon laser) 
provides a beam of monochromatic light 48 which is 
expanded and collimated by lenses 50 (e.g., spatial filter 
and collimating lens from Newport Corp.) to form a 
collimated beam 51 which is passed through the optical 
modulator 45. The light beam 51 is diffracted by the 
image of the Young’s fringes on the modulator 45, with 
the partially diffracted output beam 52 from the modu- 
lator being passed through a polarizer 53 (e.g., sheet 
polarizer, Oriel Corp.) and lens 55 (e.g., Nikkor Macro 
105) to produce an image at the focal plane 56 which 
corresponds to the optical Fourier transform of the 
Young’s fringes image on the modulator 45. Another 
lens 57 (e.g., microscope objective, such as 3 mm focal 
length achromat from Newport Corp.) positioned at the 
image plane 56 expands the light pattern at the plane 56 
and images it onto a video camera 60 (e.g., RCA vid- 
icon without lens) which produces a coaxial cable out- 
put signal on an output line 61 corresponding to a raster 
scan of the image striking the video camera. The result- 
ing analog signal is converted to digital data in an ana- 
log to digital converter 62 (which may be within the 
frame grabber) and provided to the microcomputer 
system 40 for further processing. 

The processing of each frame from the video camera 
60 consists of DC subtraction to remove the center peak 
in the image recorded on the camera 60 (as explained 
below), peak finding (preferably using a line integral 
algorithm) to find the positions of these peaks, data 
storage to store the locations of the two peaks for each 
frame, film advance to send a signal to the stepper 
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motor 31 to move the film to a new point for interroga- 
tion, and then a repeat of the steps above. 

A typical output from the microcomputer as dis- 
played on a video monitor for one laser probe beam 
position is shown in FIG. 4, illustrating a central DC 5 
spot and two off axis correlation peaks corresponding to 
the input fringe angle and displacement. A convenient 
low pass filtering technique is a frame subtraction, of a 
low frequency only convolution, from the high fre- 
quency frame modulated beam convolution. Results 10 
using this frame subtraction technique are illustrated in 
FIGS. 5-8 below. In each of the figures, the figure 
marked “A” is the input to the optical modulator 
45 — the image captured by the video camera 36—, and 
the figure marked “B” is the output of the microcom- 15 
puter 40 after low frequency subtraction. FIGS. 5-8 
show three different fringe orientations (all at the same 
spatial frequency), and their respective convolutions. 
FIG. 8 shows a lower frequency fringe spacing. Each of 
these images were produced as “phantom” fringes 20 
which were calculated numerically, rather than from a 
picture captured by the tv camera 36, and then used to 
modulate the image on the modulator 45 to allow pre- 
cise control over the orientation. The phantom fringes ^ 
allow characterization of the system for minimum re- 
solvable fringe frequencies, as well as allowing the ef- 
fects of contrast ratio and synthetic speckle noise to be 
determined. 

For further comparison purposes, the fringe pattern 
of FIG. 5A was simulated numerically and digitally 
transformed using a 256 x 256 two dimensional fast 
Fourier transform (FFT) routine. A comparison of a 
slice through the correlation peak of the digital trans- 
form with a similar slice for the processor output of the 35 
present invention (via the video camera 60 and an eight 
bit frame grabber of conventional design) is shown in 
FIG. 9. For this comparison, the digital transform was 
rescaled to match the video transform peak value and 
frequency. The two dimensional FFT took about 40 4 0 
milliseconds on a Cray XMT with the FFT algorithm 
written in assembly language and optimized for the 
Cray’s vector operation. The figure illustrates the rela- 
tive conformance of the width of the video correlation 
peak with respect to the “idealized” digital correlation 45 
peak width. 

The following is a source code listing, written in C 
and labeled t42b.c, for carrying out the image process- 
ing utilizing the subroutine library available with the 
frame grabber board (Imaging Technologies FG-100). 50 

The basic elements of the program are: 

Set-up phase, establishing values, limits, etc., from 
lines 1-49. This portion sets the step size and range of 
motion for the stepper motors, loads the reference 
image used for image subtraction, sets up the input 55 
look-up table to stretch the pixel values of the image if 
necessary, and finally implements the peak-finder to 
find the center of the reference image. 

Loops, (lines 50-56 and 61-65) set up the repetition 
scheme, to allow the evaluation of the image point by 60 
point in a raster scan type of pattern. The loops incre- 
ment the location of the point being interrogated by 
moving the stepper motors in the appropriate direction 
the appropriate amount. 

Snap (line 57) captures the current image on camera 65 
60 on to the computer frame grabber board. Since the 
picture is now in memory, changes to the image on the 
camera have no effect until another snap is issued. 
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Frame-subtract (line 58) calls a procedure (lines 
69-89) which subtracts a reference image from the re- 
cently captured image. The image subtraction allows 
filtering the low frequency content out of the captured 
image. 

Movex (line 59) calls a procedure (lines 312-330) 
which drives the x-axis stepper motor to move the film 
a certain amount in a given direction. 

Peak-finder (line 60) calls a procedure (lines 91-310) 
which finds the location of the autocorrelation peaks in 
the image. The peak finder uses an integral sum search 
method in which line integrals of the pixel values are 
calculated at 20 pixel intervals in both x and y directions 
throughout the image. The maximum integral in the x 
direction and the maximum integral in the y direction 
are assumed to give coordinates close to the higher of 
the two expected autocorrelation peaks. A second 
search is initiated based on this approximate location, 
using line integrals in x and y directions at 1 pixel inte- 
grals. To save calculation time, the lines are truncated 
to 60 pixels long. The maximum x and y integrals from 
this search are considered the exact location of the 
autocorrelation peak. In order to find the second peak, 
the first peak is blocked out of the image, and the entire 
search procedure is reinitiated. The locations of the two 
peaks are stored on a ramdisk for speed. 

After the peaks are found, (note the film has already 
been moved) a new image is snapped, and the process is 
repeated for each interrogation point in the field. 

These are the basic steps in the program. Alterna- 
tively, the DC subtraction can be done optically with a 
high pass filter instead of by image subtraction. 

It is understood that the invention is not limited to the 
particular embodiments described herein, but encom- 
passes all such modified forms thereof as come within 
the scope of the following claims. 

What is claimed is: 

1. A method of image processing for use in measuring 
the velocity of particles flowing in a fluid, wherein 
images of the particles have been recorded on optical 
medium in an optical double exposure showing the 
displacement of the particles over a time increment, 
comprising the steps of: 

(a) projecting light through a portion of the medium 
onto a screen to form an image of Young’s fringes: 

(b) modulating a second beam of light in a two dimen- 
sional array in accordance with the image formed 
on the screen to produce a diffraction pattern in the 
light so modulated; 

(c) forming an image of the diffraction pattern. 

2. The method of claim 1 wherein steps (a), (b) and (c) 
are repeated in a sequence with the light in step (a) 
projected through a different portion of the medium for 
each repetition to cover the entire image recorded on 
the medium. 

3. The method of claim 1 further including the step of 
providing an electrical signal corresponding to the 
image of the diffraction pattern. 

4. The method of claim 2 further including the step of 
providing an electrical signal corresponding to the 
image of the diffraction pattern at each of the portions 
of the medium through which light is projected. 

5. The method of claim 3 wherein the step of provid- 
ing an electrical signal the signal corresponds to the 
light intensity of the diffraction pattern in a raster scan 
of the image. 

6. The method of claim 5 further including the steps 
of converting the electrical signal to digital data and 
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analyzing the data to find the frequency and orientation 
of the diffraction pattern at the portion of the medium 
through which the light was projected. 

7. The method of claim 1 wherein the step of project- 
ing light through a portion of the medium is carried out 5 
by activating a laser to direct a narrow beam of mono- 
chromatic laser light through the portion of the me- 
dium, and wherein the step of modulating a second 
beam of light, the second beam of light is provided by 
activating a laser to provide a collimated beam of mono- 10 
chromatic light which is modulated in the two dimen- 
sional array in accordance with the image formed on 
the screen. 

8. The method of claim 1 wherein the step of forming 
an image of the diffraction pattern is carried out by 15 
focusing the modulated second beam of light onto a 
video camera. 

9. The method of claim 1 wherein the step of modu- 
lating a second beam of light is carried out by passing a 
beam of monochromatic light through a two dimen- 20 
sional spatial light modulator having a two dimensional 
matrix of cells of variable transmissiveness, each of the 
cells being controlled to attenuate the light impinging 
thereon in relation to the light intensity of the corre- 
sponding part of the image on the screen which mani- 25 
fests the Young’s fringes. 

10. Apparatus for image processing for use in measur- 
ing the velocity of particles flowing in a fluid, wherein 
images of the particles have been recorded on an optical 
medium in an optical double exposure showing the 30 
displacement of particles over a time increment, com- 
prising: 

(a) means for projecting light through a portion of the 

optical medium onto a screen to form an image 
manifesting Young’ fringes; 35 

(b) means for modulating a second beam of light in a 
two dimensional array in accordance with the 
image formed on the screen to produce a diffrac- 
tion pattern in the light so modulated; 

(c) means for viewing the diffraction pattern and 40 
forming an electrical signal corresponding thereto. 

11. The apparatus of claim 10 further including means 
for analyzing the electrical signal corresponding to the 
diffraction pattern to determine the magnitude and ori- 
entation of the fringes at the portion of the recorded 45 
medium through which the light was projected. 

12. The apparatus of claim 10 further including means 
for holding the optical medium and for moving the 
optical medium in a pattern in step wise motion so that 
light is projected through the medium over the entire 50 
area of the medium constituting the recorded informa- 
tion. 

13. The apparatus of claim 10 wherein the means for 
modulating includes a video camera mounted to record 
the image formed on the screen and to provide an elec- 55 
trical signal corresponding thereto and a spatial light 
modulator connected to receive the signal from the 
video camera and having a two dimensional array of 
cells which are controlled to pass light therethrough in 
relation to the intensity of a portion of the image formed 60 
on the screen. 

14. The apparatus of claim 13 wherein the spatial 
light modulator comprises a rectangular array of liquid 
crystal cells. 

15. The apparatus of claim 10 wherein the means for 65 
projecting light through a portion of the medium in- 
cludes a laser which projects a narrow beam of mono- 
chromatic light which is directed to pass through a 
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small portion of the medium and with the light pro- 
jected therefrom striking the screen, and wherein the 
screen comprises a translucent screen through which 
the image formed on the screen can be viewed. 

16. The apparatus of claim 13 wherein the means for 
modulating a second beam of light further includes a 
laser providing an output beam of monochromatic light 
as the second beam of light and means for expanding 
and collimating the output beam from the laser and 
passing the collimated beam through the spatial light 
modulator. 

17. The apparatus of claim 14 including a panel of 
polarizer material mounted to receive the light modu- 
lated by the spatial light modulator and to polarize the 
light, and further including lenses to receive light 
passed through the polarizer material and to focus the 
light to form an image of the diffraction pattern. 

18. Apparatus for image processing for use in measur- 
ing the velocity of particles flowing in a fluid, wherein 
images of the particles have been recorded on an optical 
medium in an optical double exposure showing the 
displacement of particles over a time increment, com- 
prising: 

(a) means for projecting light through a portion of the 
optical medium onto a screen to form an image 
manifesting Young’fringes; 

(b) means for modulating a second beam of light in a 
two dimensional array in accordance with the 
image formed on the screen to produce a diffrac- 
tion pattern in the light so modulated; 

(c) means for viewing the diffraction pattern and 
forming an electrical signal corresponding thereto; 
and 

(d) means for analyzing the electrical signal corre- 
sponding to the diffraction pattern to determine the 
magnitude and orientation of the fringes at the 
portion of the recorded medium through which the 
light was projected. 

19. THe apparatus of claim 18 further including 
means for holding the optical medium and for moving 
the optical medium in a pattern in step wise motion so 
that light is projected through the medium over the 
entire area of the medium constituting the recorded 
information. 

20. The apparatus of claim 18 wherein the means for 
modulating includes a video camera mounted to record 
the image formed on the screen and to provide an elec- 
trical signal corresponding thereto and a spatial light 
modulator connected to receive the signal from the 
video camera and having a two dimensional array of 
cells which are controlled to pass light therethrough in 
relation to the intensity of a portion of the image formed 
on the screen. 

21. The apparatus of claim 20 wherein the spatial 
light modulator comprises a rectangular array of liquid 
crystal cells. 

22. The apparatus of claim 18 wherein the means for 
projecting light through a portion of the medium in- 
cludes a laser which produces a narrow beam of mono- 
chromatic light which is directed to pass through a 
small portion of the medium and with the light pro- 
jected therefrom striking the screen, and wherein the 
screen comprises a translucent screen through which 
the image formed on the screen can be viewed. 

23. The apparatus of claim 20 wherein the means for 
modulating a second beam of light further includes a 
laser providing an output beam of monochromatic light 
as the second beam of light and means for expanding 
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and collimating the output beam from the laser and 
passing the collimated beam through the spatial light 
modulator. 

24. The apparatus of claim 21 including a panel of 
polarizer material mounted to receive the light modu- 5 


10 

lated by the spatial light modulator and to polarize the 
light, and further including lenses to receive the light 
passed through the polarizer material and to focus the 
light to form an image of the diffraction pattern. 

***** 


10 


15 


20 


25 


30 


35 


40 


45 


50 


55 


60 


65 



UNITED STATES PATENT AND TRADEMARK OFFICE 

CERTIFICATE OF CORRECTION 

Page 1 of i3 

5,011,278 

Apri l 30, 1991 

Patrick V. Farrell 
It is certified that error appears in the above-identified patent and that said Letters Patent is hereby 

corrected as shown below: 


PATENT NO. 

DATED 

INVENTORY) 


Column 1, line 32, "Pulse" should be --Pulsed--. 

Column 5, line 50, after "FG-100)." add --FILE LISTING 
for t42 b . c 


1 

2 

3 

4 

5 

6 

7 

8 
9 

10 
1 1 
1 2 

13 

14 
1 5 
16 

17 

18 
1 9 
20 
21 
22 

23 

24 

25 

26 


#include 
^include 
# i nc lude 
U i n c l u d e 
tt i nc l ude 


< s t L i b . h> 
<stdio.h> 
<dos . h> 
"itexIOO.h" 
"stdtyp.h" 


peak_finder() ; 
f rame_subt ractO; 
movexO ; 
movey ( ) ; 
int nx,ny; 


void 
void 
void 
void 
long 
int dirx,diry; 

FILE * f open ( ) ,* f p; 
m a i n ( ) 

■C 

unsigned long delta_ms(); 

unsigned int n t ot x , nt ot y , i x , i y , L l o, l h i ; 

int err; 

float s t epx , snx , s t epy , sny ; 

char name[20],commentC200],contC10] ; 

sethdw(Ox300,OxAOOUG,2,1) ; 

setdim(1024,512,12); 

initialized ; 

sclear(0); 

setreg(0x10,0) ; 

printfl" File name of reference image file: 


\ n " ) 
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27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 

41 

42 

43 

44 

45 

46 

47 

48 

49 

50 

51 

52 

53 

54 

55 

56 

57 

58 

59 


scant ( ”%s" ,name) ; 

err=readim(512, 0,512, 512, name, comment); 
printfC'Step size in x direction? \n"; 
scant C"%f", 8 stepx) ; 

printfC'Step size in y direction? \n"); 
scanfl"%f", 8 stepy) ; 

printfC "Number of steps in the x direction? \n") 
scanfC"%d", 8 ntotx) ; 

printfC "Number of steps in the y direction? \n") 
scanfC"%d",Sntoty> ; 

printfC "Lut stretch lo and high values? \n"); 

scanf C "%d %d " ,8 l lo, 8 L h i ) ; 

snx=stepx/5.08e-5; 

nx=Clong)snx; 

sny=stepy/5.08e-5; 

ny= C long) sny; 

fp=fopenC"f:pfiv.dat","w" ); 
snap C 1 ) ; 
peak_finderC) ; 

lut_stretchCINPUT,0,llo,Lhi); 
dir x = 0 ; 
d i ry= 0 ; 
iy = 0 ; 

while C iy < ntoty ) -C 
if C dirx > 0 ) 
d i r x = 0 ; 
else 
di rx = 1 ; 
i x = 0 ; 

while C ix < ntotx ) -C 
snapC 1 ) ; 

frame_subtractC) ; 
movexCdi rx) ; 
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60 

61 

62 

63 

64 

65 

66 

67 

68 

69 

70 

71 

72 

73 

74 

75 

76 

77 

78 

79 

80 
81 
82 

83 

84 

85 

86 

87 

88 

89 

90 


peak_f i nde r ( ) ; 

+ + i x ; 

> 

movey ( d i ry ) ; 

+ + iy; 

> 

p r i nt f ( " \ n D on ' t forget to copy data file from 
ramdisk (f:)! "); 

> 

void frame subtractO 
-C 

int i , j , a n s, e r r , mag, ang l e, of l ag ; 

static unsigned px a [ 5 1 2] , px bC 5 1 2 ] , p x c [ 5 1 2] ; 

i=0; 

while ( i < 480 ) -C 

rhline(0,i,512,pxa) ; 
rhline(512,i,512,pxb) ; 
j=0; 

while (j < 512) -C 

if( pxbUj] > pxaCj] ) 
px c [ j 3 = 0 ; 
else 

pxcCj3=pxa[j3-pxb[j3; 

+ + j ; 

> 

w h l i ne ( 0 , i , 5 1 2 , px c ) ; 

+ +i ; 

> 
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91 void 

peak finderO 

92 

•c 

93 

unsigned int pxsum,pysum,ixmax / iymax; 

94 

int i,j,err,pval / ict / irx r iry,ii,ictx,icty; 

95 

int ictxmax,ictymax,iimin,iimax,ijmin,ijmax; 

96 

int iclear,jclear,peakx1,peaky1 ; 

97 

static unsigned pxa[51 23 ,px[51 23 ,pyL51 23 ; 

98 

i c t = 0 ; 

99 

i r x = 1 5 ; 

100 

while (irx < 496) { 

101 

i c t= i c t + 1 ; 

102 

rvlineCi rx,15,450,pxa) ; 

103 

ii = 1; 

104 

px sum=0 ; 

105 

while (ii < 450) -C 

106 

pxsum=pxsum+pxaCii3 ; 

107 

+ + i i ; 

108 

> 

109 

pxCi ct3=pxsum; 

110 

i rx= i rx+20; 

1 1 1 

> 

112 

i c t x= i c t ; 

113 

i c t = 0 ; 

114 

i ry = 1 5 ; 

1 1 5 

while ( iry < 464 ) { 

116 

i c t= i c t + 1 ; 

1 1 7 

rhline(15,iry,482,pxa); 

1 1 8 

ii = 1; 

1 1 9 

py sum = 0 ; 

1 20 

while ( ii < 482 ) -C 

121 

pysum=pysum+pxaCi i3 ; 

122 

+ + i i ; 
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123 

> 

124 

py C i ct] = py sum; 

125 

i ry= i ry + 20 ; 

126 

> 

127 

i c t y= i c t ; 

128 

i i = 2 ; 

129 

i x m a x = 0 ; 

130 

while ( ii < ictx ) { 

131 

if ( px [ i i 3 > ixmax 

132 

i xma x = p x C i i ] ; 

133 

ictxmax=i i ; 

1 34 

> 

1 35 

+ + i i ; 

136 

> 

137 

i i = 2 ; 

138 

i y ma x = 0 ; 

139 

while C ii < icty ) { 

140 

if ( pyCii] > iymax 

141 

i ymax=py [ i i D ; 

142 

i c t y ma x= i i ; 

143 

> 

144 

+ + i i ; 

145 

> 

1 46 

printf( M fi rst 

1 47 

i imax=( i ctxmax-1 )*20+20; 

1 43 

iimin=(ictxmax-1)*20-20; 

149 

i j m a x = ( i ctymax-1 )*20+20; 

1 50 

i jmi n=( i ctymax-1 )*20-20; 

151 

i c t = 0 ; 

152 

i = i i m i n ; 

153 

while (i <= iimax) { 


) -C 


) -C 


thru %d %d ", i ctxmax, i ctymax) 
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1 54 

ict=ict+1; 

155 

rvline(i,ijmin,40,pxa) ; 

156 

i i =0; 

157 

px sum=0 ; 

158 

while (ii < 40) -C 

1 59 

pxsum=pxsum+pxaCii]; 

160 

+ + i i ; 

161 

> 

162 

px L i ctD=pxsum; 

163 

+ +i ; 

164 

> 

165 

ictx=ict; 

166 

i=i jmi n; 

167 

i c t = 0; 

168 

while ( i <= ijmax ) -C 

169 

i ct= i ct + 1 ; 

170 

r h L i ne ( i imin,i,40,pxa) 

171 

i i = 0 ; 

172 

py sum=0 ; 

173 

while ( ii < 40 ) -C 

174 

pysum=pysum+pxaCi 

175 

+ + i i ; 

176 

> 

1 77 

pyC i ct]=pysum; 

1 78 

+ +i ; 

179 

> 

180 

i c ty= i c t ; 

181 

i i = i ; 

182 

i xmax=0 ; 

1 83 

while ( ii < ictx ) -C 

184 

if ( pxCii] > ixmax ) 
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185 

i xmax = px[ 

186 

i ct xma x = C 

187 

> 

1 88 

+ + i i ; 

1 89 

> 

1 90 

i y ma x = 0 ; 

191 

ii = i; 

192 

while ( ii < icty ) 

193 

if ( py C i i 3 > 

1 94 

iymax=py[ 

195 

i ctymax= ( 

1 96 

> 

1 97 

+ + i i ; 

198 

> 

1 99 

peakx 1 =i ctxmax; 

200 

peaky 1 =i ctymax; 

201 

iclear=ictxmax-40; 

202 

j c lear=i ctymax-40; 

203 

aclear(iclear,jclea 

204 

i c t = 0 ; 

205 

i r x = 1 5 ; 

206 

while (irx < 496) { 

207 

i c t = i c t + 1 ; 

208 

rvline(irx,15. 

209 

ii = i; 

210 

px sum = 0 ; 

21 1 

while (ii < 45 

212 

pxsum=pxs 

213 

+ + i i ; 

214 

> 

215 

p x [ i c 1 3 = px sum; 
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216 

i r x= i r x + 2 (J ; 

2 17 

> 

21 a 

ictx=ict; 

219 

i ct = 0 ; 

2 2u 

i ry = 1 5 ; 

221 

while ( iry < 464 ) { 

222 

i ct = i ct + 1 ; 

226 

rhline(15,iry,482,pxa) ; 

2 24 

i i = 1 ; 

225 

py sum=0 ; 

226 

while ( ii < 482 ) { 

227 

pysum=pysum+pxaCi i ] ; 

228 

+ + i i ; 

229 

> 

2 3D 

py C i ct]=pysum; 

231 

i r y = i r y + 2 0 ; 

232 

> 

233 

i c t y= i c t ; 

234 

i i=2; 

235 

i x m a x = U ; 

236 

while ( ii < ictx ) { 

237 

if ( pxCii] > ixmax ) -C 

238 

i xmax=px[ i i ] ; 

23y 

i ct xma x= i i ; 

240 

> 

241 

+ + i i ; 

242 

> 

243 

ii=2; 

244 

i y ma x = 0 ; 

245 

while ( ii < icty ) -C 

246 

if ( pyCii] > iymax ) { 
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247 

iymax=py[i i 3 ; 

248 

i c tymax= i i ; 

249 

> 

2 5 0 

+ + i i ; 

251 

> 

252 

i imax=( i ctxmax-1 )*20+20; 

253 

i i mi n= ( i ctxmax-1 ) *20-20; 

254 

i j max=( i ctymax-1 )*20 + 20; 

255 

i jmi n=( i ctymax-1 )*20-20; 

256 

i ct = 0 ; 

257 

i = i i m i n ; 

258 

while (i <= iimax) { 

259 

i ct = i ct + 1 ; 

260 

rvline( i jmi n,40,pxa) ; 

261 

i i = 0 ; 

262 

px sum = 0 ; 

263 

while (ii < 40) { 

264 

pxsum=pxsum + pxa[i ill; 

265 

+ + i i ; 

266 

> 

267 

pxEi ct]=pxsum; 

268 

+ + i ; 

269 

> 

270 

ictx=ict; 

271 

i = i j m i n ; 

272 

i ct = 0 ; 

273 

while ( i <= ijmax ) { 

274 

i ct=i ct + 1 ; 

275 

rhlineli imin,i,40,pxa) ; 

276 

i i =0 ; 

277 

py sum = 0 ; 
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278 

while ( ii < 40 ) -C 

279 

pysum=pysum+pxa[i i ] ; 

280 

+ + i i ; 

281 

> 

282 

pyCi ct]=pysum; 

283 

+ +i ; 

284 

> 

285 

i c t y= i c t ; 

286 

i i =1 ; 

287 

i xmax=Q; 

288 

while ( ii < ictx ) -C 

289 

if ( p x C i i □ > ixmax ) { 

290 

i xma x = px [ i i ] ; 

291 

ictxmax=i i - 1 + i imin; 

292 

> 

293 

+ + i i ; 

294 

> 

295 

i y ma x = G ; 

296 

i i = 1 ; 

297 

while ( ii < icty ) -C 

298 

if ( pyCii] > iymax ) { 

299 

i ymax=py[ i i 3 ; 

300 

ictymax=i i-1+i jmin; 

301 

> 

302 

+ + i i ; 

303 

> 

304 

f pr i nt f ( f p, " %d %d %d %d \ n" , pea kxl , pea kyl , i c t xma x, 
i ctymax) ; 

305 

printf("%d %d %d %d \ n" , pea kxl , peakyl , i c t x ma x, i c t y ma x ) 

306 

hlclear(0,peaky1 ,512,255) ; 

307 

vlclearlpeakxl ,0,480,255) ; 

308 

hlclear(0,ictymax,512,255) ; 
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309 v l c l ea r ( i c t x ma x , 0 ,480 , 2 5 5 ) ; 

310 > 


void 


void 


signed int a d d r e s s , x x x , y y y , x p , x n ; 
ig int i ; 

Ires s = 81 6 ; 

<= add re s s + 1 5 ; 

add re s s + 1 2 ; 

: 4* d i rx+1 ; 

= 4 * d i rx ; 
tp(xxx,128); 
tp(yyy,xn> ; 

); 

i l e ( i < nx ) { 

outp(yyy,xp) ; 
outp(yyy,xn) ; 

+ + i ; 

> 

tp(yyy,48) ; 

/ey ( ) 

signed int a d d r e s s , x x x , y y y , y p , y n ; 

ig int i ; 

iress=816; 

<=address+1 5 ; 

^=address+1 2 ; 

=8*d i ry+2 ; 

- 8 * d i r y ; 



UNITED STATES PATENT AND TRADEMARK OFFICE 

CERTIFICATE OF CORRECTION 


PATENT NO. 

DATED 

INVENTOR(S) 


5,01 1 ,278 
Apri L 30, 1 991 
Patrick V . Farrell 


Page 12 of 13 


It is certified that error appears in the above-identified patent and that said Letters Patent is hereby 
corrected as shown below: 

340 outp(xxx,128); 

341 outp(yyy,yn); 

342 i =0 ; 

343 w h i l e ( i < ny ) { 

344 outp(yyy,yp) ; 

345 outp(yyy,yn) ; 

346 ++i; 

347 > 

348 outp(yyy,48); 

349 > 


Column 6, line 45, "fringes:" should be --fringes;- 
Column 7, line 35, " Y oung ' f r i n g e s , " should be 
--Young's fringes;-- 

Column 7, line 67, "projects" should be —produces- 
Column 8, line 26, "Young'fringes," should be 
--Young's fringes,--. 

Column 8, line 39, "THe" should be --The--. 
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